Anomalous radial expansion in central heavy-ion reactions 
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The expansion velocity profile in central heavy-ion reactions in the Fermi energy domain is examined. 
The radial expansion is non-Hubblean and in the surface region it scales proportional to a higher 
exponent [a > 1) of the radius. The anomalous expansion velocity profile is accompanied by a 
power law nucleon density profile in the surface region. Both these features of central heavy-ion 
reactions disappear at higher energies, and the system follows a uniform Hubble expansion (a ~ 1). 



I. INTRODUCTION 

In heavy-ion collisions, the stage of compression and 
heating is followed by the expansion of nuclear matter. 
Expansion dynamics as a collective motion of excited 
matter is characterized by certain space-momentum cor- 
relations and has been well ascertained by experimental 
data such as the collective fiow. In particular, in central 
heavy-ion collisions with the beam energy ranging from 
the Fermi energy to almost 200 GeV/A, the radial flow 
is clearly manifested through the flattening of the trans- 
verse spectra with the particle mass and this effect, as 
expected, is stronger for heavier systems [0. 

The collective expansion scenario is important also in 
other issues. In studying the quantum statistical correla- 
tions which describe the space-time characteristics of ex- 
panding systems, one can infer the information about the 
freeze-out configuration. The recent finding here is that 
the size parameters of an effective source are determined 
not only by the geometrical length scale which measures 
the region of homogeneity ||] but also by the thermal 
length scale which is related to the region in the coor- 
dinate space from which identical particles with similar 
momenta may emerge Alongside with the tempera- 
ture and the freeze-out time, the thermal size of a source 
is determined by the velocity gradients in this source. 
The thermal length dominates the correlation function if 
the geometrical length scale is sufficiently large Q . Sim- 
ilarly, the question about the particle velocity profile at 
the freeze-out configuration arises when one considers the 
formation of light nuclei in the framework of a coalescence 
model |]. 

In almost all papers devoted to this subject, the veloc- 
ity profile is assumed to be linear : 
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where v and f are three- or two-dimensional vectors for 
the spherically or cylindrically symmetric expansion, re- 
spectively. The relation (|l|) is prompted by an analytical 



scaling solution of the equations of non-relativistic hy- 
drodynamics with the ideal gas equation of state for a 
slowly expanding fireball 1^,^. In fact, Eq. (|l|) is a con- 
sequence of a regular motion governed by the continuity 
equation when a fluid does not influence the expansion 
rate. This relation is well known in cosmology where the 
Hubble constant (7|J| : R{t)jR{t) « 65 km/s/Mpc, at 
the right-hand side of Eq. (Q), characterizes the expan- 
sion of homogeneous and isotropic galaxies § . 

It is important to note that Eq. (|^) need not be valid 
in general. As was shown by Dumitru while the lon- 
gitudinal expansion of the flreball three-volume is inde- 
pendent of the energy density of the fluid, the transverse 
collective motion in the case of (3-|-l)-dimensional ex- 
pansion may couple the expansion rate to the properties 
of the fluid, i.e., to the equation of state. In particu- 
lar, the hydrodynamic solution for a flreball expanding 
in the longitudinal and transverse direction with a possi- 
ble flrst order hadronization phase transition affects the 
three- volume expansion rate on the hadronization hyper- 
surface H]. 

Analyzing the experimental data on nuclear multifrag- 
mentation within an extended statistical microcanonical 
model which takes into account an interplay of the radial 
expansion with a non-spherical shape of the fragmenting 
nucleus |]l0|, it turned out to be necessary to postulate 
a non-Hubblean velocity proflle v{r) for fragments in the 
freeze-out volume : 



v{r) = Wo 



r 



(2) 



with a in between 1.5 and 2 pO| , to explain the exper- 
imental charge-number dependence of the mean kinetic 
energy of fragments in central Xe + Sn collisions at 50 
MeV/A. The answer to the question : why the exponent 
a of the radial expansion differs from 1 (a = 1 yields the 
Hubble expansion), cannot be given within the statistical 
multifragmentation model. An application of the hydro- 
dynamics for this kind of problem is also questionable 
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because the nucleon density at the freeze-out configura- 
tion is low and, moreover, the dynamical processes at 
short time scales are not correctly described. For that 
reason, in this work we study both the particle veloc- 
ity profile and the particle density profile in central HI 
collisions using the framework of a nonlocal quantum ki- 
netic theory To gain an insight into the dynamics 
of collective expansion of small fermionic systems such 
as the atomic nuclei, the kinetic approach which uses the 
quasiparticle interaction as input and takes into account 
consistently the two-particle correlations is probably 
more reliable, even though the dynamical formation of 
clusters is absent in this approach. 

The paper is organized as follows. In Sect. II. A, 
the main ingredients of the nonlocal quantum kinetic 
approach are presented. The time evolution of central 
Ta+Au collisions at 33 MeV/A and 60 MeV/A is studied 
in Sect. II. B by looking at the transversal and longitu- 
dinal profiles of the nucleon velocity, the nucleon density 
and the proton to neutron ratio. The expansion veloc- 
ity profile is discussed in more details in Sect. II. C , 
separately for bulk and surface particles. The qualitative 
evolution of the radial expansion profile with the collision 
energy is compared in Sect. II. D with the dynamical tra- 
jectories of excited system in the temperature - particle 
density plane. The possible consequences of the long- 
range tail in the particle density on the small momenta 
behavior of the Bose-Einstein correlations in discussed in 
Sect. II. E. Finally, Sect. Ill summarizes main results of 
the paper. 



II. THE KINETIC APPROACH 

A. The nonlocal quantum kinetic equation 

The observables of interest are : the particle density 
n{r,t), the current density J{r,t) and the kinetic energy 
density E{r, t), which can be expressed by the one - par- 
ticle phase - space distribution function /(p, r, t) as fol- 
lows : 



n(r, t) 
J{r,t) 
E{r,t) 



dp 



(27r)3 

dp 
(2^ 



3 P f{p,r,t) 



(3) 
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The one - particle distribution function obeys a nonlocal 
Boltzmann - Uehling - Uhlenbeck (BUU) kinetic equation 

m : 



dfi ^ 9£i 9/i _ 9£l 9/l _ 

dt dk dr dr dk ^ 
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with the Enskog-type shifts of the arguments 

fi = f{k,r,t) 

f2 = f{p,r~A2,t) 

f3 = f{k-q-AK,r- A3, t- At) 

f4 = f{p+q-AK,r-A4,t-At) . 



(5) 



/3/4(l-/l)(l-/2)-(l-./3)(l-/4)/l/2 , (4) 



The arguments of the effective scattering measure V are 
centered in all A - shifts. The quasiparticle energy e con- 
tains the mean field as well as the correlated self energy. 
The shifts or displacements are a compact form of gra- 
dient corrections and ensure that the conservation laws 
contain both the mean-field and the two-particle correla- 
tions. In particular, the momentum and the energy gain 
arises from the finite duration of collisions . All shifts 
in are proportional to derivatives of the scattering 
phase shift ||ll|,|l^,|l^ and have been calculated for realis- 
tic nuclear potentials [|l5|. When neglecting these shifts, 
one recovers the usual BUU scenario. 

It should be noted that using the nonlocal BUU kinetic 
equations for the description of the proton spectra 
in central Xe + Sn collisions at 50 MeV/A leads to a 
significant enhancement of the high energy tail [p^ and 
a better agreement with the experimental data than ob- 
tained using the standard BUU equations. 



B. The evolution plots 

The result of the nonlocal BUU scenario for the reac- 
tion Ta + Au at 33 MeV/A can be seen in Figs. |^ and 
^. Let us concentrate first on the corresponding velocity 
and density profiles (the first and the second column in 
Fig. |) and the arrows characterizing the mass momen- 
tum (the first column in Fig. ^ . One sees that at around 
40 fm/c the nuclei start to squeeze out the matter side- 
wards (the first column in Fig. 2), what is characterized 
by the momentum focusing at both sides perpendicular 
to the beam direction , predominantly in peripheral re- 
gions. The surface matter is stopped and bounced back 
in longitudinal direction during the times 20 — 80 fm/c 
(see the first column in Fig. |^) . The inner (bulk) matter 
exhibiting a quite clear spatial boundary (see the second 
column in Fig. |^) is still moving inwards. This leads at 
short time scales (~ 40 fm/c) to an enhancement of mat- 
ter density. The strong velocity gradient, which is seen 
at 60 fm/c, disappears at about 100 fm/c and the inner 
matter comes to a rest. The recoil of the splashing matter 
at the surface and the attractive mean field force start 
effectively to reaccelerate an inner matter towards the 
center of mass. Since the surface particles are still accel- 
erating in the outward direction, therefore there is a zone 
of matter in between which comes to a rest. This evolu- 
tion leads to the dumb-bell structure in the transversal 
density profile at 120 - 240 fm/c. 

The different behavior of the surface matter and the 
bulk matter leads to the development of a nonlinear ve- 
locity profile in the surface region, which can be seen in 
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the log - log representation of the angular averaged ve- 
locity in the second column of Fig. ||. For t > 80 fm/c, 
the velocity - radius scaling with an asymptotically sta- 
ble coefficient asurf — 1.75±0.05 (see the second column 
in Fig. 1^) appears definitely in the surface region. In this 
region, the particle density drops nearly as a power law : 
n(r) ~ r~^, with the asymptotically stable (for t > 120 
fm/c) coefficient : /3 ~ 3 ± 0.2. As can be seen from the 
N/Z - ratios in the third column in Fig. the Coulomb 
interaction expels protons from the surface in the early 
stage of the evolution. In particular, the proton rms ra- 
dius is larger than the neutron one. This effect becomes 
weaker at later times {t ~ 200 fm/c) but, nevertheless, 
it survives indicating that proton and neutron distribu- 
tions are different in the surface region. This may lead 
to different source temperatures and temperature gradi- 
ents for protons and neutrons and, hence, to different 
quantum statistical corrections for protons and neutrons 
in the interferometry experiments. This possibility has 
been suggested in the phenomenological analysis of the 
asymmetric reaction Ar + Au at 30 MeV/A [ [17[ . 

These two unusual effects : the nonlinear velocity pro- 
file with asurf G [1-5,2.0], and an approximately power 
law fall-off of the particle density with /3 ~ 3, charac- 
terize the transitional region in the central collisions of 
symmetric HI collisions in the Fermi energy domain. 

For later times {t > 200 fm/c), one sees the forma- 
tion of an oblate configuration which is connected to the 
inversion of the velocity of the inner matter and to the ac- 
cumulation of the density. In agreement with earlier ob- 
servations [|l8|,|l9| , we see that the formed hot and nearly 
fused matter is not spherically symmetric. We would 
like to remark that this deformation is specific for ener- 
gies around the Fermi energy and, moreover, is impact- 
parameter dependent. At nonzero impact parameters, 
the shape of the matter distribution becomes prolate due 
to the spectator matter keeping its initial direction of mo- 
tion and also due to the angular momentum effects. 

The evolution picture shown in Figs. |l| and ^ is es- 
sentially changed at higher bombarding energies. At 
Eiab/A = 60 MeV (see Figs. | and |), the time inter- 
val where the Coulomb force counter-balance the nuclear 
forces is becoming very short {e.g., see the plots for t = 40 
fm/c) and the system enters very fast in the phase of a 
smooth radial expansion. In this case, the velocity - ra- 
dius scaling can be well approximated asymptotically by 
a single exponent abuik = Q-surf — 1- The particle den- 
sity is expanding almost uniformly and no characteristic 
power law dependence is seen in the surface region. The 
proton excess in the surface region is less pronounced 
than at Eiab/A = 33 MeV/A (see Fig. 1) and Z/N ~ 1 
at later times. One expects that with increasing bom- 
barding energy in symmetric HI reactions the effective 
source parameters for protons and neutrons become close 
to each other. 



C. The expansion velocity 

Let us now discuss the dependence of the expansion ve- 
locity on the radius in more details. As was noted above, 
two different slopes can be distinguished at lower colli- 
sion energies. We call the 'inner' and the 'outer' parts 
of the density profile the 'bulk' and the 'surface' regions, 
respectively. They are separated here at around R = 10 
fm. Therefore, we plot in Fig. ^ the time dependence 
of the exponent a in these two regions for different bom- 
barding energies. 

Let us start with the reaction at Eiab/A ~ 33 MeV. 
The very first stage of the collision, where the nuclei are 
slightly overlapping, is characterized by similar values of 
the exponents a in the bulk and the surface regions. This 
feature continues until the surface particles begin to be 
evaporated. At around this time, the surface develops 
a much steeper velocity gradient with asurf as large as 
^2. At the same time, the bulk matter velocity profile 
is quite smooth and even abuik changes the sign [ po| . 
After this overshooting of the surface exponent during a 
time interval 100 — 200 fm/c, the bulk matter develops 
the radial expansion with the coefficient abuik which is 
approaching asurf ■ Note, that even at this stage the mass 
current is characterized by the nonlinear scaling. Hence, 
using the Hubble ansatz for the radial flow in the analysis 
of experimental data in the Fermi energy domain is not 
justified. 

At lower energies {Eiat/A = 15 MeV in Fig. 5), one ob- 
serves that the difference between radial flow patterns in 
the surface and bulk regions becomes even stronger than 
at Eiab/A = 33 MeV. Consequently, higher values of the 
exponent asurf are reached asymptotically. Moreover, 
we see that a giant resonance with a period of T = 60 
fm/c, or an energy of 211 /T = 20.6 MeV, is excited. This 
can be considered as a complete fusion event. 

If one proceeds to energies higher than the Fermi en- 
ergy, one sees that the general trend of evolution is con- 
served but the deviation between surface and bulk matter 
expansion patterns decreases, i.e., the coefficients abuik 
and asurf become close one to another. One can see also 
that the maximum value of the exponents a is reached 
faster than at lower energies, and the exponents a are 
becoming close to 1 at late times . The nuclear system 
at this high excitation energy is expanding continuously. 
Due to the higher initial velocity, there is practically 
no inversion of the velocity profile and no time-periodic 
structures like giant resonances are excited. The Hub- 
blean expansion pattern is reached faster and no anoma- 
lous behavior [usurf > 1) is observed for energies higher 
than 90 MeV/A. 

D. Dynamical trajectories 

To elucidate the connection of the anomalous velocity 
profile with the multifragmentation, let us characterize 



3 



an instantaneous state of the system in terms of the av- 
erage nucleon density n and the temperature T. In order 
to define a global time dependent temperature T{t) we 
adopt the Fermi liquid relation pi| : 



E{t) = -Ep{t) 



EcoW {t) 



iEpit) 



nty 



(6) 



where the global kinetic energy E{t), the collec- 
tive energy ii'coii and the Fermi energy Ep{n) = 
(37r^n/2)^/'^/2m are given by a spatial integration of the 
local quantities (3) : 



E{t) 



Epit) = 



Eco\\{i) — 



drE{r, t) 
dr n{r, t) 

dr Ep {n{r,t)) n{r,t) 



dr 



(7) 



dr n{r, t) 



mn{r, t) 



dr n{r, t) 



It is more problematical to define a density. We will 
present here the two possibilities. The first one is to 
consider the density of matter inside the evolving mean 
square radius. The other possibility is to consider the 
density of matter contained in the static separation of 
bulk matter by the radius R found earlier in the velocity 
profile. 

The dynamical trajectories in the {T ~ n) - plane for 
different collision energies are shown in Fig. 6 for the 
bulk matter (i? < 10 fm) . Let us first look at the static 
density definition inside the bulk region. For Eiab/A = 
15 and 33 MeV, the system evolves inside the spinodal 
region and, therefore, is mechanically unstable. In the 
time interval between 150 fm/c and 200 fm/c, the sys- 
tem is in a configuration with n « 7io/3 and T « 7 — 8 
MeV, which are the typical values for the nuclear multi- 
fragmentation. At higher energies, the freeze-out density 
of an evolving system is shifted towards lower densities 
and finally it ends in a gaseous phase. If one chooses 
to look at the interior part of the system within the ra- 
dial size of the mean squared radius (dotted lines in Fig. 
6), the main difference is seen in the initial stage of the 
evolution where larger nucleon densities are reached but, 
at the same time, they are passed through much faster. 
The freeze-out configurations are practically the same as 
when the above static definition of bulk matter is consid- 
ered. 



E. The long tail of the density distribution 

Let us now discuss the particle density profile in more 
details. As noted in Fig. 2, the decrease of the parti- 
cle density in the surface region is algebraic (power law) 
rather than exponential. In Fig. 7, we plot the time de- 
pendence of the exponent (3 which is extracted from the 
power law fit : n{r) ~ r~^, for i? > 10 fm. 

There are here two distinct behaviors. For Eiab/A =15 
and 33 MeV, after an initial build-up of the surface re- 
gion as characterized by the decrease of (3 in time, the 
value of exponent /? stabilizes asymptotically. The limit- 
ing value of /3, called f3urm decreases with bombarding en- 
ergy from I3iim ~ 3.5 at Eiab/A = 15 MeV to (3iim ~ 3.1 
at Eiab/A = 33 MeV. Actually, the most interesting 
long-range region of the nucleon density ends at about 
Eiab I A ~ 50 MeV and the smallest value reached is about 
= 2. In the energy interval 15 MeV < Eiab/A < 50 
MeV, the anomalous power law tail of the nucleon den- 
sity accompanies the anomalous profile of the expansion 
velocity in the surface region, and both effects are lead- 
ing to the non-Gaussian shape of the emitting source and 
the anomalous short-range correlations. 

For energies Eiab/A > 60 MeV in the Hubblean ex- 
pansion regime, the exponent (3 decreases monotonously 
in time and no asymptotically stable surface region with 
the power law dependence is seen. At later times , one 
sees however the appearance of a new expansion regime 
corresponding to the negative values of the exponent (3. 
This indicates the formation of a shell-like structure in 
the system. The formation time of the shell-like struc- 
ture decreases rapidly with increasing bombarding en- 
ergy {tsh ^ 280 fm/c at Eiab/A = 60 MeV and t,h ^ 170 
fm/c at Eiab/A = 90 MeV). The /3 values found are not 
strongly negative, indicating that the expanding shell is a 
very diffused object due to the long mean-free path of nu- 
cleons in the kinetic approach. Similar unusual solution 
of spherically expanding scaling hydrodynamic has been 
used in the analysis of the Bose-Einstein correlations . 

The specific shape fluctuations in the long-range region 
can be a source of the power-law Bose-Einstein correla- 
tions. This problem has been discussed by Bialas [^3|j24[ 
in the context of strong interaction physics at relativis- 
tic energies. The long tail of the particle density for the 
systems produced in the symmetric heavy-ion collisions 
in the Fermi energy domain, may result in the unusual 
quantum statistical correlations. The TV-particle inter- 
ferometry reduced densities can be written as : 



Dpfiki, . . . ,A:jv) 



= dxi . .. dxN ^ exp[i(a;ipai 

V per 



+ . . . + XNPaM)] X 



X npfixi, . . . ,xn) 



Nl / dxni{x) 



N 



(8) 



where is the A^-point density of emitting sources and 
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the sum runs over all permutations of the indices ai. The 
above formula supposes the incoherent emission from the 
source and neglects the final state interaction. In the case 
of uncorrelated emission in space-time : 

nNixi, . . . ,xn) = ni{xi) ■ ■ -niixN) , 

the A'^-particle cumulant can be written using the Fourier 
transform of the source density ni{x) |23] : 



UN 



(fci, . . . , /cat) = ^ ni{ki - kaj - ■ ■ ni{kN - ka^) , (9) 



where the sum runs over all permutations of the indices 
ai with ai ^ i. If the source density has a power law tail 

mix) :^x-f^ -x^-'^ , (10) 

its Fourier transform shows also a power law in some 
range of small momenta [E3| : 



ni{k) 



(11) 



Thus there is the relation between the power law tail in 
the source density distribution and the power law in the 
two-particle Bose-Einstein correlations which are given 
in terms of the Fourier transform of the source density : 



C2{q) ^ \q\ 



-27 



(12) 



Similarly, the higher order cumulants Ci{q) {i = 3, . . .) 
are expected also to have a power law dependence on the 
rescaling of momenta with an index 17. 

A particularly interesting case of the Bose-Einstein cor- 
relations, increasing as the power of \q\, corresponds to 
7 > (/3 < D) in Eq. (p^). In the case of central 
heavy-ion collisions, our analysis suggests that this effect 
could be seen in symmetric systems for energies : 35 MeV 
< Eiab/A < 50 MeV. 



III. THE SUMMARY AND OUTLOOK 

The dynamical behavior of heavy-ion reactions in the 
Fermi energy domain during the first 200 fm/c is clearly 
associated with a nonlinear (a > 1) radial velocity pro- 
file (2). The existence of such a non-Hubblean radial 
expansion at the freeze-out configuration was postulated 
by Le Fevre et al. [|oj in the analysis of experimental 
kinetic energies of fragments in the Xe + Sn reaction 
at Eiab/A — 50 MeV in the framework of the statisti- 
cal microcanonical model. The present studies using the 
nonlocal kinetic theory show that indeed such an unusual 
radial flow velocity profile is a plausible freeze-out con- 
figuration for symmetric heavy-ion collisions. One can 
roughly determine the instant of time when the compres- 
sion turns into the expansion by looking at the crossing 
point of the bulk ahuik and surface asurf scaling expo- 
nents. While at the compression stage these two ex- 
ponents are almost equal, the expansion stage exhibits 



larger exponent for the surface matter than for the bulk 
matter and asurf takes values significantly larger than 
1. This dynamical behavior of the surface matter dis- 
appears for energies significantly higher than the Fermi 
energy. 

For central heavy-ion collisions in the Fermi energy 
domain, we see that the expansion stage is character- 
ized by very small values of the exponent a in the bulk 
{ctbuik & [0, 1]), indicating its slow evolution and possi- 
ble mechanical instabilities. At these energies, the sys- 
tem spends a long time in the spinodal region, what may 
result in its multifragmentation decay. With increasing 
collision energy, the system passes quickly through this 
unstable region or, at even higher energies, its thermody- 
namic trajectories go above the spinodal region. In this 
case, the multifragmentation due to the spinodal insta- 
bilities is hardly possible . For details see ^l[ 

The nonlinear radial velocity profile in the surface re- 
gion is accompanied by the long tail of the particle den- 
sity. Both these effects may have important consequences 
on the quantum statistical correlations and their evolu- 
tion with the bombarding energy. The effect of the ra- 
dial expansion on the two-particle correlations have been 
studied assuming the linear scaling solution (a = 1) of 
the scaling hydrodynamics [pyi7|]. It was found that the 
expansion makes the effective radius of the two-particle 
correlation functions smaller than the geometrical size of 
the source. We expect this effect to be even stronger in 
the presence of the nonlinear radial expansion flow, lead- 
ing to even stronger discrepancy between the effective 
radius and the geometrical radius. Moreover, the com- 
monly used Gaussian approximation for the source shape 
is certainly hazardous in this non-Hubblean expansion 
regime. The existence of the power law tail in the parti- 
cle density can in turn lead to the power law two-particle 
Bose-Einstein correlations at small relative momenta of 
the particles in the range of collisions energies 35 MeV 
< Eiab/A < 50 MeV. It should be stressed that the ratio 
Z/N in the surface region is strongly different from 1 , in 
particular at early collision times and at low bombarding 
energies {Eiab/A < 60 MeV). This specific effect in the 
Fermi energy domain, which leads to an effective increase 
of the proton rms radius and to the change in the tem- 
perature gradients in proton/neutron source, may have 
measurable consequences in the quantum statistical cor- 
relations for protons and neutrons. This effect has been 
studied by Helgesson et al. |17| assuming a linear scaling 
solution (a = 1) of hydrodynamics for asymmetric HI 
reactions at 30 MeV/ A. Our results suggest that the si- 
multaneous description of n and p spectra, as well as nn 
and pp correlation functions may require different source 
parametrizations for neutrons and protons, though the 
detailed dynamics for HI reactions at Eiab/A ~ 30 MeV 
is very different from those assumed by Helgesson et al. 
pj{ . It is interesting to notice , that the nonlocal kinetic 
theory predicts the appearance of the solution somewhat 
similar to the linear scaling solution (a = 1) of hydrody- 
namics with /? < [p2[ for higher energies bombarding 
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energies {Eiat/A > 60 MeV. 

The different behavior of the surface and bulk matter 
before equihbration can also be of importance for the de- 
scription of super nova where a surface - like ring of mat- 
ter (crust) is expanding with enormous velocities and is 
clearly separated from the remaining bulk matter collaps- 
ing back into neutron stars. The clearest experimental 
observation of the expanding shell-like structure noted 
above comes from stellar astronomy. The envelope ma- 
terial ejected by the stars forms an expanding shell of 
gas that is known as a planetary nebula. The space-time 
evolution of these objects is in many aspects similar to 
the considered evolutions : a < 2, P > and a ~ 1, 
P < 0. The latter solution can be successfully simulated 
by a scaling solution of the non-relativistic hydrodynam- 
ics H. 

We suggest that the anomalies found in the kinetic ex- 
pansion reflect the nature of effective interactions among 
the elementary constituents of the system. Their man- 
ifestation is twofold. Firstly, the interplay between the 
repulsive Coulomb interaction and the attractive mean 
field results in the formation of a rather sharp surface 
of the system. Secondly, the evolution of the bulk mat- 
ter is not a simple uniform expansion of a homogeneous 
ideal fluid. In the unstable spinodal region, the local 
interaction of quasiparticles in the bulk phase affects no- 
ticeably the subsequent evolution of the system. In this 
respect, one should stress again an attractive possibil- 
ity of implementing both the nonlinear velocity profile 
and the algebraic long-range density tail into the anal- 
ysis of the Bose-Einstein correlations. These effects are 
not only important for the nuclear multifragmentation 
process but also for the hadron interferometry at ultrarel- 
ativistic collisions where the deconfinement phase transi- 
tion can have a strong influence on the expansion stage. 
This is of a particular interest for the statistical mixed 
phase equation of state which gives the crossover type of 
the deconfinement phase transition and allows for a small 
admixture of unbound quarks at the freeze-out point . 

An interesting example which illustrates an important 
difference between the linear and quadratic scalings of 
the expansion velocity with the radius can be found in 
cosmology 0. If there is an attractive force decreasing 
as an inverse power of the radius, the equation of motion 
for a radially symmetric matter follows from the total 
energy h which reads : 

Assuming a homogeneous matter density n, the mass is : 
m — AttuR^/S. In the case of the escaping matter {h — 
0), one gets : 

I 3-5 

what corresponds to taking : 



in Eq. (^. We see that 6—1 for the Coulomb or gravi- 
tational forces and, hence, the Hubble expansion (a = 1) 
follows in these cases. However if a = 2, as found in the 
surface region of nuclei formed in central HI collisions at 
around the Fermi energy (see Fig. 5), the above rela- 
tions lead to a string- like force with 6 — ~1. Since this 
force is used to describe the confinement in the effective 
theories motivated by the Quantum Chromodynamics, 
its occurrence as a consequence of the dynamical behav- 
ior is worth of attention. As was shown above in the 
solution of nonlocal kinetic equations, the interplay be- 
tween Coulomb and mean field leads to such a string-like 
behavior for surface particles. This is here a clear non- 
equilibrium effect. The above discussed oscillation in the 
time-dependence of the a exponent may be considered 
as a possible manifestation of this effective dynamical 
string-like force in the surface region. 

Perhaps the best way to demonstrate the existence of 
both the non-Hubblcan radial expansion and the alge- 
braic long-range tail of the particle density, would cor- 
respond to finding a non-Gaussian deformation of the 
source in the Bose-Einstein interference experiments for 
central, symmetric heavy-ion reactions in the narrow 
range of collisions energies (35 MeV < Eiab/A < 50 
MeV). In the same narrow range unusual scale - depen- 
dence of the many-particle correlations at small momenta 
should be induced via the same Bose-Einstein quantum 
interference effect if the source has an algebraic (power 
law) density tail (Sect. II. E). On the other hand, as dis- 
cussed in pffl , the kinematical observables related mainly 
to the intermediate mass fragments, provide an indepen- 
dent and sensitive measure of the velocity profile in the 
deformed, expanding source. To put together all these 
different pieces of evidence into the circumstantial proof 
for the non-Hubblean collective expansion of nuclear ag- 
gregates, remains a difficult and exciting challenge for the 
future experimental and theoretical studies of heavy-ion 
collision in the Fermi energy domain. 
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FIG. 1. Tlie time evolution of central To -I- Au collisions at Eiat/A = 33 MeV in the nonlocal kinetic model The first 
column represents the velocity profile. Both transversal (the a;-direction) (the dashed line) and the longitudinal (the z-direction) 
(the solid line) profiles are shown. The second column shows the longitudinal (the solid line) and the transversal (the dashed 
line) density profiles in units of fm~^. The third column presents the ratio of proton and neutron densities normalized to the 
initial Z/A - ratio for both the longitudinal (the solid line) and the transversal (the dashed line) projections. 
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FIG. 2. The time evolution of central Ta + Au collisions at Eiab/A = 33 MeV in the nonlocal kinetic model [|l^. Plots in 
the first column show the {x — z) - density cut. The mass momenta are shown by arrows. The second column shows the log - 
log plot of the angular averaged modulus of the expansion velocity versus the radius. The slope of the straight line fit of the 
surface matter expansion profile for > 10 fm is indicated at each plot (see also Fig. ^). The third column shows the log - 
log plot of the angular averaged nucleon density versus the radius. The slope of the straight line fit of the surface profile is 
indicated at each plot as well. 
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FIG. 3. The same as in Fig. 1 but for Eiab/A = 60 MeV. 
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FIG. 5. The exponent a of the power law velocity profile (2) with respect to the radius v oc r" for different lab energies. The 
dotted lines show the surface matter behavior and the solid lines depict the bulk matter behavior. 
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FIG. 6. The time evolution of the temperature versus the density for the sphere of a constant radius {R < 10 fm) (the sohd 
hne) corresponding to the distinction between the bulk and surface regions in Fig. 1, and for the sphere given by the mean 
square radius of the evolving system (the dotted line). The thick solid lines give the limit of the spinodal region for an infinite 
nuclear matter. 
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